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ABSTRACT 

We computed proper motions of a selected sample of globular clusters projected on 
the central bulge, employing CCD images gathered along the last 25 years at the ESO- 
NTT, ESO-Danish and HST telescopes. We presented a method to derive their proper 
motions, and a set of coordinate transformations to obtain 3D Galactic velocity vectors 
of the clusters. We analysed 10 globular clusters, namely Terzan 1, Terzan 2, Terzan 4, 
Terzan 9, NGC 6522, NGC 6558, NGC 6540, AL 3, ES0456-SC38 and Palomar 6. For 
comparison purposes we also studied the outer bulge cluster NGC 6652. We discuss the 
general properties of the proper-motion-cleaned Colour-Magnitude Diagrams, derived 
for the first time for most of them. A general conclusion is that the inner bulge 
globular clusters have clearly lower transverse motions (and spatial velocities) than 
halo clusters, and appear to be trapped in the bulge bar. 

Key words: Galaxy: bulge - Globular clusters: individual: NGC 6522, NGC 6558, 
AL 3, HP 1, ES0456-SC38 (Djorgovski 2), NGC 6540 (Djorgovski 3), Terzan 1, Terzan 
2, Terzan 4, Terzan 9, Palomar 6, NGC 6652 


1 INTRODUCTION 

Globular clusters (GCs) in the Galactic bulge preserve in 
their spatial distribution and orbital evolution essential in¬ 
formation to probe the early the early formation stages of 
the Galaxy central parts. The combination of dynamical 
properties of globular clusters, with their ages and chemical 
composition, provides a new tool to investigate the bulge 
stellar populations, and to build a consistent scenario of the 
Galactic bulge formation. 

Stars and globular clusters in the Galactic halo present 
very elliptical orbits, with low angular momentum, while 
disk objects show circular orbits with high angular momen¬ 
tum, and small vertical velocity. The bulge instead shows an 
intermediate angular momentum, with higher perpendicular 
velocities than disc stars, but not going as deep into the halo 
as genuine halo stars. For the field bulge stars, kinematics 
reveals two different behaviours (Babusiaux et al. 2010): (1) 
a metal poor component, enriched in [Mg/Fe], consistent 
with an isotropic rotating population belonging to an old 
spheroid, and (2) a metal rich one with a vertex deviation 
consistent with that expected from a population with orbits 
supporting a bar. The spatial separation of these two com¬ 
ponents is not well determined. More recently, Babusiaux et 
al. (2014) and Zoccali et al. (2014) carried out large surveys 


on stars in the Galactic bar, inner disk and central bulge, by 
measuring radial velocities and metallicities from GaT lines, 
and traced the bar with more precision. Of interest to us is 
the conclusion that there is a concentration of metal-poor 
stars, that could have formed before the bar instability, and 
could have been trapped by the bar when it formed, staying 
confined in the innermost regions of the boxy bulge (Babu¬ 
siaux et al. 2014; Di Matteo et al. 2014). Their conclusion 
is that the exact distribution of metal-poor stars, and their 
connection with the formation history of the thick disk, in¬ 
ner halo and the bar need further investigations. 

Most of previous efforts in bulge cluster proper motion 
measurements were carried out using HST data as those e.g. 
by Zoccali et al. (2001) for NGC 6553, Feltzing & John¬ 
son (2002) for NGC 6528, Bellini et al. (2013) for NGC 
6338 and NGC 6441, and Kuijken & Rich (2002) for field 
stars. Ground-based large telescopes equipped with Multi- 
Conjugate Adaptive Optics also allowed proper-motion- 
cleaning to be performed, such as those carried out by Fer¬ 
raro et al. (2009) for Terzan 5, and Ortolan! et al. (2011) for 
HP 1. 

Calculations of Galactic orbits, based on proper motions 
and radial velocities, for some inner Galaxy clusters were 
carried out by Dinescu et al. (2003). 

In previous works we studied the inner bulge globular 
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clusters, by means of Colour-Magnitude Diagrams (CMDs) 
and spectroscopy (e.g. Ortolani et al. 1995; Barbuy et al. 
1998, 2009). We now have an unprecedented archive of CCD 
images taken since the 90s with European Southern Obser¬ 
vatory (ESO) NTT and Danish telescopes, and the Hubble 
Space Telescope (HST). Through the combination of this 
archive data with more recent CCD images, this data set 
allows cluster proper motions to be derived, and therefore 
proper-motion-cleaned CMDs and kinematics (e.g. Ortolani 
et al. 2011). 

We focus on moderately metal-poor GCs ([Fe/H]~ —1) 
with blue horizontal branch, projected at low galactic lat¬ 
itude, because they might be the oldest population in the 
Galaxy (e.g. Barbuy et al. 2009). They might belong to the 
same generation of stars as the central RR Lyrae (Dekany et 
al. 2013), the latter identified as an old and spheroidal com¬ 
ponent of the bulge. We also study a few more metal-rich 
GCs. 

Our main goal is to derive accurate (~0.5 mas/yr) ab¬ 
solute proper motions. Kinematic properties of the sample 
clusters should allow bulge members and halo intruders to be 
distinguished. The proper-motion results, coupled with ra¬ 
dial velocities, allow 3D-orbit determinations in the Galactic 
potential (e.g. Ortolani et al. 2011). The orbits will be used 
to constrain different stellar-population components. In a 
forthcoming paper the orbits will be presented in detail. 

In this work we derive the initial state vector for a total 
of ten GCs together with a recompilation of their metallic- 
ities. Sect. 2 describes the selection of the sample. In Sect. 
3 we report the observations. In Sect. 4 proper motions are 
derived. In Sect. 5 we present an analysis of errors, the re¬ 
sults are presented in Sect. 6 and conclusions are drawn in 
Sect. 7. 


2 SAMPLE SELECTION 

2.1 Metallicity and spatial distribution 

The metallicity distribution of the Milky Way globular clus¬ 
ters is shown in Figure [T] where metallicities are adopted 
from Bica et al. (2006), in some cases updated with more 
recent spectroscopic analyses (Barbuy et al. 2009, 2014). 

Figure [T] shows a Gaussian distribution centred around 
[Fe/H] ~ —1.6, representing the metal-poor population, and 
a Gaussian distribution around the value [Fe/H] ~ —0.55, 
corresponding to the metal-rich population (e.g. Morgan 
1959; Kinman 1959; Zinn 1985; Cote 1999). A moderately 
metal-poor sub-population characterized by [Fe/H] ~ —1.0 
is evident (see also Cote 1999). 

With this criterion in mind we subdivided the globular 
clusters into three sub-samples: 

• metal-poor GCs (95 objects): [Fe/H] < —1.1 

• metal-rich GCs (42 objects): [Fe/H] > —0.9 

• moderately metal-poor GCs (13 objects): —1.1 ^ 
[Fe/H]- ^ -0.9. 

In this work we prioritize clusters with a blue horizontal 
branch together with a metallicity of [Fe/H]~ —1.0, includ¬ 
ing NGC 6522, NGC 6558, AL 3, ES0456-SC38, NGC 6540. 
We also study two mainly red horizontal branch and mod¬ 
erately metal-poor clusters, Terzan 1 and Palomar 6, two 



Figure 1. Metallicity distribution of the Milky Way globular 
clusters, updated with recent literature. Dashed lines are gaussian 
fits for the metal-rich and metal-poor populations. 


metal-poor globular clusters projected on the inner bulge 
(Terzan 4, Terzan 9), and one metal-rich (Terzan 2). For 
comparison purposes we also included the outer bulge NGC 
6652, with proper motion values taken from Sohn et al. 
(2014). 

The log of observations is presented in Table l??l In 
Table 2 we report the cluster distances from the literature. 
This Table shows a good agreement among authors, which 
gives confidence on the distance values. 


2.2 The distance of the Galactic centre 

The distance of the Sun from the Galactic centre (R©) is a 
crucial parameter that we use to establish a cluster initial 
state vector. 

One of the most common methods adopted to determine 
the value of this parameter is the determination of the mean 
distance of Galactic globular clusters (Bica et al. 2006), or 
the distribution of RR Lyrae stars (Dekany et al. 2013). 
Another method is based on the distribution of field clump 
stars (Bobylev et al. 2014; Nishiyama et al. 2006; Cao et al. 
2013). Finally, Gillessen et al. (2009), and Do et al. (2013) re¬ 
ported recent values from the orbit of stars around SGR A*, 
and derived a geometrical distance. A recent review of these 
methods and results was presented by Malkin (2013). The 
distances of the Sun to the Galactic center essentially vary 
from 7.5 kpc to 8.5 kpc. In the present paper we considered 
both 7.5 kpc, for the reason that this distance is consistent 
with the distribution of bulge globular clusters, and 8.3 kpc, 
as reported by Brunthaler et al. (2011), Schonrich (2012), 
Reid et al. (2014) and Hon & Han (2014), among others. In 
Table |3] we give the proper motion values, and Galactic dis¬ 
tance components and velocity vectors, computed for both 
Rq = 7.5 kpc, and Rq = 8.3 kpc (see Sect. 4). 


2.3 Distribution of clusters in the Galaxy 

Figure [2] shows the Galactocentric-distance distribution of 
the three metallicity sub-samples of GCs. Figures [3] and 
[4] show the Galactic distribution of the three sub-samples. 
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Table 1. Log of observations of first and second epochs. 


Cluster 

Telescope 

Instrument 

Date 

Filter 

Exposure 

Seeing 

Terzan 1 1st 

Danish 

0.47”/pix 

06/1990 

I 

180s 

1.1 


Danish 

0.47”/pix 

06/1990 

Gunn z 

120s 

1.1 

Terzan 1 2nd 

NTT 

EFOSC2 0.24”/pix 

05/2012 

I 

10s 

1.1 


NTT 

EFOSC2 0.24’7pix 

05/2012 

Gunn z 

180s 

1.1 

Terzan 2 1st 

NTT 

SUSI 0.13”/pix 

05/1994 

V 

30,720s 

0.8, 0.9 


NTT 

SUSI 0.13”/pix 

05/1994 

I 

60,420s 

0.8, 0.8 

Terzan 2 2nd 

NTT 

EFOSC2 0.24”/pix 

06/2011 

V 

60,300,20s 

1.1,1.1,1.1 


NTT 

EFOSC2 0.24”/pix 

06/2011 

I 

40,180,10s 

1.1,1.1,1.1 

Terzan 4 1st 

NTT 

SUSI 0.13”/pix 

05/1994 

V 

60,600s 

0.55, 0.55 


NTT 

SUSI 0.13”/pix 

05/1994 

I 

60,300s 

0.4,0.4 


NTT 

SUSI 0.13”/pix 

05/1994 

Gunn z 

120s 

0.5 

Terzan 4 2nd 

NTT 

EFOSC2 0.24”/pix 

05/2012 

I 

30s 

0.6 


NTT 

EFOSC2 0.24”/pix 

05/2012 

V 

30s 

0.6 

Terzan 9 1st 

Danish 

DFOSC 0.39”/pix 

07/1998 

V 

60,900,60s 

1.3,1.3,1.3 


Danish 

EFOSC2 0.39”^ix 

07/1998 

I 

60,90s 

1.1,1.1 

Terzan 9 2nd 

NTT 

DFOSC 0.24”/pix 

05/2012 

V 

30,90s 

0.5,0.7 


NTT 

EFOSC2 0.24”/pix 

05/2012 

I 

10,20s 

0.55,0.6 

NGC 6522 1st 

HST 

WFPC2 

09/1995 

F439W 

50,160,160s 

— 


HST 

WFPC2 

09/1995 

F555W 

10,50s 

— 

NGC 6522 1st 

Danish 

0.47”/pix 

06/1992 

V 

60,480s 

1.2,1.3 


Danish 

0.47”/pix 

06/1992 

I 

20,300s 

1.2,1.2 

NGC 6522 2nd 

NTT 

EFOSC2 0.24”/pix 

05/2012 

V 

30,60s 

0.55,0.75 


NTT 

EFOSC2 0.24”/pix 

05/2012 

I 

15,30s 

0.5,0.7 

NGC 6558 1st 

NTT 

EMMI 0.35”/pix 

06/1993 

I 

7,120s 

1.2,1.2 


NTT 

EMMI 0.35”/pi^ 

06/1993 

V 

10,180s 

1.2,1.2 

NGC 6558 2nd 

NTT 

EFOSC2 0.24’/pix 

05/2012 

V 

90,300s 

0.6,0.6 


NTT 

EFOSC2 0.24Vpix 

05/2012 

I 

60,180,20s 

0.5,0.55,0.55-0.6 

NGC 6540 1st 

Danish 

0.47”/pix 

06/1990 

V 

60s 

1.3 


Danish 

0.47”/pix 

06/1990 

I 

15s 

1.0 

NGC 6540 2nd 

NTT 

EFOSC2 0.24”/pix 

05/2012 

V 

30s 

0.9 


NTT 

EFOSC2 0.24”/pix 

05/2012 

I 

10s 

0.9 

AL 3 1st 

Danish 

DFOSC 0.39”/pix 

03/2000 

V 

60,180s 

1.2,1.2 


Danish 

DFOSC 0.39”/pix 

03/2000 

I 

10,40s 

1.2,1.2 

AL 3 2nd 

NTT 

EFOSC2 0.24”^ix 

05/2012 

V 

30,60,240,60s 

0.6 for all 


NTT 

EFOSC2 0.24”/pix 

05/2012 

I 

10,30,180s 

0.6 

ES0456-SC38 1st 

NTT 

SUSI 0.13”/pix 

05/1994 

I 

60s 

0.9 


NTT 

SUSI 0.13”/pix 

05/1994 

V 

180s 

0.8 


NTT 

SUSI 0.13”/pix 

05/1994 

Gunn z 

60s 

0.8 

ES0456-SC38 2nd 

NTT 

EFOSC2 0.24”/pix 

05/2012 

V 

900,60s 

1.0,1.0 


NTT 

EFOSC2 0.24”/pix 

05/2012 

I 

30s 

0.8 

Palomar 6 1st 

NTT 

EMMI 0.35”/pix 

06/1993 

V 

900s 

1.4 


NTT 

EMMI 0.35”/pix 

06/1993 

I 

300s 

1.4 

Palomar 6 2nd 

NTT 

EFOSC2/HR 0.12”/pix 

05/2013 

Gunn z 

300s 

0.9 


NTT 

EFOSC2/HR 0.12”/pix 

05/2013 

I 

300s 

1.1 


The coordinates of the clusters are expressed in the inertial 
right-handed Galactocentric frame of reference, in which the 
a;-axis points towards the Sun. The model of the Galactic 
bar and of the spiral pattern have been chosen according to 
Pichardo et al. (2004) and Pichardo et al. (2003), respec¬ 
tively. In Figs. [2l [3] and [4] we used a Galactocentric distance 
of Rq = 7.5 kpc. 

We identified the following features of the three sub¬ 
samples. Metal poor clusters are distributed in a nearly 
spherical volume with a radius of about 30 kpc, consis¬ 


tently with a halo population. The spatial distribution of 
metal-rich clusters is more centrally concentrated than the 
metal-poor ones, as expected. The moderately metal-poor 
sub-sample is confined within the inner 5 kpc of the Milky 
Way. 

In order to identify the clusters located in the Galactic 
bulge, we selected a region on the celestial sphere included 
between |Z| ^ 18° and |6| ^ 7.5°. The metallicity distribution 
of the 39 clusters located in this region of the sky is shown in 
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Table 2. Distances from the literature (all values in kpc). 


Cluster 

Bica+06 

Valenti+07 

Harris96 

Barbuy+98 

other 

Reference 

Terzan 1 

6.2 

6.6 

6.7 

4.90 

5.2 

Ortolan! et al. (1999a) 

Terzan 2 

8.7 

7.4 

7.5 

6.64 

7.7 

Ortolan! et al. (1997) 

Terzan 4 

9.1 

6.7 

7.2 

7.28 

8.0 

Ortolan! et al. (2007) 

Terzan 9 

7.7 

5.6 

7.1 

— 

4.9 

Ortolan! et al. (1999b) 

NGG 6522 

7.8 

7.4 

7.7 

6.05 

6.2 

Barbuy et al. (1994) 

NGG 6558 

7.4 

— 

7.4 

- 

6.3 

Rich et al. (1998) 

NGG 6540 

3.7 

5.2 

5.3 

3.02 

3.5 

B!ca et al. (1994) 

AL 3 

— 

— 

6.5 

— 

6.0 

Ortolan! et al. (2006) 

ES0456-SG38 

6.7 

7.0 

6.3 

3.02 

5.5 

Ortolan! et al. (1997) 

Palomar 6 

7.3 

— 

5.8 

- 

8.9 

Ortolan! et al. (1995) 

NGG 6652 

9.6 


10.0 

— 

9.3 

Ortolan! et al. (1994) 


22 
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Figure 2. Galactocentric distance (in kpc) distribution of the 
three sub—samples of globular clusters. 

Figure [S] The moderately metal-poor clusters show a peak 
comparable to that of the metal-rich clusters. 

In Figure [5] we show the projection of their coordinates 
on the Galactic plane and on the x-z plane, where a distance 
to the Galactic center of Rq = 7.5 kpc was adopted. They 
appear centrally concentrated, and possibly trapped within 
the bar. 


3 OBSERVATIONS 

The first epoch archival data are based on the ESO NTT 
telescope, obtained in 1993 and 1994, the ESO Danish 1.5m 
telescope in 1990, 1998 and 2000, and the Hubble Space Tele¬ 
scope (HST). The latter observations were acquired in 1995, 
under the proposal GO 9065, PI: S. Djorgowski, and the 
photometry is available as reported in Piotto et al. (2002(0. 

The NTT 1993 data used the EMMI spectrograph/focal 
reducer equipped with the GCD Loral ESO #34, of 

^ http://www.astro.unipd.it/globulars/ 



Figure 3. Spatial distribution of globular clusters projected on 
the Galactic plane on the x—y plane, viewed from the North Galac¬ 
tic Pole. 


2048x2048, giving a O'.'35/pixel scale. The NTT 1994 im¬ 
ages employed the SUSI camera equipped with the GGD 
ESO #25, with 0(G3/pixel scale. The Danish 1.5m 1990 im¬ 
ages employed the direct focal camera with GGD ESO # 5, 
of 512x320 30 /rm pixels, of scale 0(#7/pixel. The Danish 
1.5m 1998 and 2000 data come from the DFOSG spectro¬ 
graph/focal reducer equipped with the Loral/Lesser G1W7 
2052x2052 15 /rm pixels, with scale O/39/pixel. 

The second-epoch observations were obtained with the 
ESO NTT-bEFOSG2 instrument, in V, I and Gunn z. The 
GGD ESO # 40 UV-flooded, of 2048x2048, with 15 /rm 
pixels was used, giving a projected scale of 0/12/pixel. The 
total held has 4(1 x 4(1. In the 2011 and 2012 observations 
the GGD was read in two by two pixel binned mode, which 
gives the best duty cycle observing time. In 2013 the full res¬ 
olution was used, because of its higher astrometric accuracy 
and higher dynamical range. 

The log of observations for the two epochs (Table l??ll 
summarize the essential information on the data. 
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Figure 4. Spatial distribution of globular clusters projected on 
the X— z plane for a selection of |Z| ^ 18° and |b| ^ 7.5°®’’°. 


Figure 6. Spatial distribution of Galactic bulge globular clusters 
overlapped on the central bar model. 



Figure 5. Metallicity distribution of the GCs projected on the 
Galactic bulge. 

4 PROPER MOTION DERIVATION 

We describe below the method applied to derive the proper 
motions of the sample. The basic idea is to transform the 
coordinates of clnster and field stars in the two different 
epochs, into a common frame of reference (main-frame) and 
compute the proper motion in terms of coordinate changes 
in the main-frame. The configuration of stars in the two 
epochs in principle should be different because of the dif¬ 
ferent proper motion of single stars. Therefore a suitable 
choice of the stars to be used to determine the coordinate 
transformation is required. We assumed that cluster mem¬ 
bers have similar proper motion and, because of the grav¬ 
itational bound, they maintain a similar configuration in 
the two epochs. We used the cluster members to determine 
the coordinate transformation in the main-frame and sub¬ 
sequently we applied the same transformation to all stars. 


both field and cluster. In this work we chose second epoch 
images as main-frame. Subtracting the value of the coordi¬ 
nates of stars in the second epoch image from those in the 
first, we expect to find (if the cluster has a significant proper 
motion) two distinct aggregations of stars in the (Aa;, Ay) 
plane: one centred around (A*, Aj/) = (0,0) (the cluster) and 
the other around another point different from the origin (the 
field). 

4.1 The technique for computing proper motions 

The raw data have been preprocessed, using MIDAS com¬ 
mands, trimming the images and subtracting the bias. Many 
tests have been carried out on the fiatfields obtained during 
the same observing runs, both inside the dome and on the 
sky twilight, and the latter proved to be better distributed. 
An average of 10 flats have been used for the flatfield correc¬ 
tion. We used the code DADPHQT (Stetson 1987) for the PSF 
photometry choosing a uniform number of bright and iso¬ 
lated stars for the average reference PSF used for fitting. A 
two degree polynomial was used across the field accounting 
for the PSF variations in the second epoch NTT-1-EFOSC2 
data. Some crowded fields needed a local subtraction of faint, 
close-by sources to the reference stars used to model the 
PSF. The fitting was performed using the ALLSTAR rou¬ 
tine which generates the relative photometry, the refined 
positions of the centroids and other parameters useful to 
clean the generated list from spurious objects or bad mea¬ 
surements. This routine creates the . als tables that we used 
as starting data for the proper motions. 

We used single-image data after several tests on the se¬ 
lected best images, as compared with the average of multiple 
observations. This is due to the peculiarity of the targets, 
that are very crowded and contaminated. The quality of the 
results is more dependent on the resolution (seeing, focus, 
tracking), and the capability to isolate the cluster CMD se¬ 
quences, rather than on the number of used frames. In ad- 
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dition, in most cases, the errors are dominated by the hrst 
epoch, where only a single colour image was available. 

The stars in common in the two images have been 
identihed and sorted out using the codes DADMATCH and 
DADMASTER. In the following analysis we use the quantities 

• (xi,j/i): coordinates of the star (in pixels) in the first 
image 

• (®2,1/2): coordinates of the star (in pixels) in the second 
image 

• (xi,t/i), (x2iJ/ 2)- coordinates of the star (in pixels) in 
the main-frame, obtained applying the transformation. 

• (Aa;, Aj/) = {x 2 — 2^1, J/2 —y\)'- difference in the position 
of the star in the main-frame in the two epochs. 


4-1.1 Field distortion and coordinate transformation 


Before proceeding further with the description of the present 
technique to compute proper motions, we propose an anal¬ 
ysis of the effects of the field distortion. In fact, the dis¬ 
tortion of star positions is a systematic effect that has to 
be treated carefully. A systematic displacement of the co¬ 
ordinates might result in a misinterpretation of the proper 
motions. A first solution to minimize this effect is to select 
stars in the central regions of the CCD, where the distortion 
is expected to be lower than in the peripheral regions. As 
discussed below, the field distortion is the reason why a con¬ 
formal transformation including shift, rotation and scale is 
not enough to transform the coordinates in the two epochs 
into the main-frame. 

We evaluated the field-distortion affecting the second- 
epoch images from the NTT. (for the first epoch we do not 
have distortion maps available). In particular, we reduced 
a set of images of a field located in the Baade’s Window 
with uniformly-distributed stars including a central frame 
and eight shifted images of the same field. The shift of the 
eight images is 60 arcsec, corresponding to approximately 
500 pixels, both in right ascension and/or in declination, 
depending on the image. Figure [7| shows the dithering pat¬ 
tern for the four images shifted both in right ascension and 
declination. For a better visualization of the dithering pat¬ 
tern we did not show the configuration of the images shifted 
only in right ascension or declination. We identified stars 
in common in all nine frames (419 stars) and computed the 
distortion vector for each star. We selected as main-frame of 
reference the system of coordinates of the central image. Ac¬ 
cording to Anderson and King (2003), in order to transform 
the coordinates of the shifted images into the main-frame 
we applied a conformal transformation on the coordinates 
of the stars in the shifted images. 

cos 9 sin 9 


x) 

t 

= X 

[yi\ 



— sin 9 cos 9 


( 1 ) 


where {xi, yi) are the coordinates of the stars in each shifted 
frame, {x\, yl) are the transformed coordinates, A is the scale 
parameter, 9 is the rotation parameter and {xo,yo) are the 
shift parameters. 

We computed the mean position of each star in the cen¬ 
tral frame of reference as follows: 
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Figure 7. Dither pattern of four of the eight images used to 
evaluate the field-distortion. 


where [xc, t/c) are the coordinates of each star in the central 
image. Our main assumption is that the mean position com¬ 
puted represents a good approximation of the distortion-free 
coordinates of the stars, justified by the fact that we ob¬ 
tained an isotropic dithering pattern. The field distortion 
has been evaluated by subtracting the coordinates in each 
shifted image to the mean coordinates. We then assumed 
that the distortion map can be well represented by a third- 
order-polynomial transformation of the form 

' Ax{x,y) = aoo + awx + aoiy + a 2 ox^ + anxy + ao 2 y^ 

asox^ + a2ix^y + ai2xy^ + aosy^ 

Ay{x, y) = 600 -b biox + boiy + 620 ®^ -b bnxy + bo 2 y^ 

, 6302 :^ + b2ix'^y + bi2xy^ -b bosy^ 

(3) 

We derived the {aij,bij) coefficients by fitting the distortion 
model to the computed distortion vectors. The top panel of 
Figure [ 8 ] shows the distortion model for NTT derived from 
this analysis, with distortion vectors magnified by a factor 
40. The central and the bottom panels show the values of 
our model of (Ax, Ay) as a function of the (x, y) coordinates 
on the CCD. The field-distortion is lower in the CCD centre 
than in the peripheral regions, where the distortion of a star 
coordinates can be up to about 2 pixels. 

In order to test the reliability of the field-distortion 
model, we applied the correction to the stellar positions 
of the 9 exposures of the Baade’s Window. We recalcu¬ 
lated the mean position of the stars in the main-frame us¬ 
ing conformal transformation and evaluated the residuals 
(dx,dt/) between the main-frame average position and those 
of the single exposures transformed on the main-frame as 
a function of the (x, y) main-frame positions (see Figure 
Ell. Note that the displacement is constant to a good ap¬ 
proximation along the whole CCD and we conclude that 
the third-order-polynomial transformation produces a good 
distortion model. Given that the clusters in our images are 
located approximately at the centre of the detector, thus 
with minimal distortions, in order to minimize the effects 
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Figure 8. Top panel:Model of the field distortion for NTT. For 
a better visualization, the distortion vectors are magnified by a 
factor 40. Central and middle panels: field distortion as function 
of the {x,y) coordinates. 
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Figure 9. Residuals of the position of the stars after correction 
for field—distortion as function of the mean coordinates. The red 
dots represent the average of the displacement. 



y (pxl) 

Figure 10. Residuals of the position of the stars in the cen¬ 
tral frame after correcting for field-distortion with a second order 
polynomial transformation as function of the mean coordinates. 


of the distortion we obtained the proper motions by select¬ 
ing field stars included within a circular area between 500 
and 1500 pixels both in x and y, centred on the cluster. 
We test if in this region of the detector the field- distortion 
can be well-represented by a polynomial transformation of 
second order. The most important reason for such a order 
decrease is that, for some sample clusters, the number of 
cluster members available to calculate coordinate transfor¬ 
mation between the two epochs is low. Decreasing the order 
will result into a more robust evaluation of the transforma¬ 
tion coefficients. Following this approach, we considered only 
the central frame and evaluated the distortion model by ap¬ 
plying the same procedure described above, decreasing the 
order of the transformation to the second power. Figure [10] 
shows the analogue results of Figure[9]for the central regions 
of the CCD after correcting the coordinates with a second 
order transformation. Also in this case the displacement is 
constant for every value of {x, y) and we can conclude that 
the distortion map in the inner regions is well represented 
by a second order polynomial model. 

Having derived a field-distortion model for all second- 
epoch exposures from NTT-I-EFOSC2, we test if a quadratic 
transformation is sufficient to transform the coordinates 
of the first-epoch images in the frame of reference of the 
second-epoch for all other instruments used in the present 
work. Explicitly: 

I x\ = ai xp -I- 02 xiyi + as yp -|- 04 xi -I- as yi + as 
\y\ = bi xP -I- b2 xiyi + bs yP -|- 64 xi + 65 yi + be 


(x2=X2 + AxM{X2,y2) 

(2/2 = 2/2 + Ai/fd(a;2,2/2) 

where {xi,yi) are the raw positions of the stars in the 
first-epoch exposure, ( 2 : 1 , i/i) are the coordinates of the 
first-epoch exposure transformed in the field-distortion- 
corrected second-epoch frame of reference, (x 2 , 2 / 2 ) are the 
raw positions of the stars in the second epoch-epoch expo¬ 
sure and (Axfd, Ayfd) are the corrections for field-distortion 
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of the second-epoch exposures. We performed the following 
tests: 

• Test 1 : Apply our model of field-distortion correction 
to the second-epoch exposure. Use quadratic transfor¬ 
mations to bring the first-epoch exposure raw positions 
on the second-epoch reference frame and measure the 
displacement between cluster members and held stars. 

• Test 2: Apply the held-distortion to both the hrst- 
epoch and second-epoch exposures. Use a conformal trans¬ 
formation to bring the hrst-epoch coordinates in the frame 
of reference of the second-epoch image and measure the dis¬ 
placement between cluster members and held stars. 

If the two methods give consistent results we should be con- 
hdent about the correctness of our approach. In the sample 
studied the only case with known held-distortion for the hrst 
epoch is NGC 6522 from HST. In particular, we adopted 
the geometric-distortion solution proposed by Anderson & 
King (2003) for WFPC2. Figure [TT] shows the results of this 
analysis. The displacement obtained from test 1 is (A®, Ay) 
= (0.24 ± 0.04,0.10 ± 0.03) pxl, while test 2 produced the 
vales {Ax, Ay) = (0.26 ± 0.05, 0.14 ± 0.04) pxl. In this case 
we evaluated the displacement by htting a double Gaussian 
to the data points. Gonsidering the errors associated to our 
measurements, we see that the results of the two methods 
are consistent and a quadratic transformation is sufhcient. 

More specihcally, we computed the transformation of 
coordinates in an iterative process. We identihed in the two- 
epoch images a circular area centred on the cluster. As dis¬ 
cussed above, the radius of such a circular area has been 
chosen in such a way as to minimize the effects of the field 
distortion. In order to estimate, to a first approximation, the 
value of the coefficients of the coordinate transformation, we 
had to identify some stars likely belonging to the cluster. We 
based our choice on some features of the CMD. In partic¬ 
ular, we assumed that most of the blue horizontal branch 
(HB) stars are cluster members. In the cases in which the 
HB was not well defined or absent, we selected the bright 
red giant branch (RGB) stars as cluster members. After the 
first fit based on the HB/RGB members, we employed an it¬ 
erative process. We applied the transformation to the whole 
sample of stars and plotted their position in the {Ax, Ay) 
plane. The cluster members were selected as the stars in¬ 
cluded within a circular area centred on the origin of the 
coordinates. The radius of the selection was chosen mak¬ 
ing sure that the selected stars were mostly located around 
the cluster evolutionary sequences of the GMD. We then 
computed again the coefficients of the transformation using 
the newly identified cluster members. We repeated this pro¬ 
cedure iteratively until the value of the coefhcients of the 
coordinate transformation converged (we verified that usu¬ 
ally 4 or 5 iterations were enough). This solution leads to 
two results: the determination of the proper motion of the 
cluster and the decontamination of the field stars from the 
cluster CMD. 


4.1.2 Gaussian interpolation 

Once the distribution of stars in the {Ax, Ay) diagram is 
known, we determined the value of proper motion (in pixels) 


Ax (pxl) 



Ax (pxl) 



Figure 11. Displacement (in pixels) of NGC 6522 derived from 
test 1 (upper panel) and test 2 (bottom panel). The black dots 
are the field stars, while the red dots are the cluster members. 


by means of an interpolation with a double Gaussian of the 
form 


f{Ax) = Aexp 
-I- Bexp 


(Aa; - Ajci)" 

2<^ci 

(Aa: - Axfieid)" 
2'’'Ield 


( 6 ) 


where {A,B), (a;ci, a:fleid) and (uci, Ufleid are the ampli¬ 
tude, the center and the dispersion of the Gaussian describ¬ 
ing the cluster stars and the field stars distribution on the 
Aa; axis respectively. The same was done for the distribution 
of stars along the Ay axis). 

Finally the differences 


AA — Aa:ci Aa:fieid 

AT = Ayd Ayfieid 


(7) 


give the cluster proper motion. 
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Figure 12. Inertial Galactocentric frame of reference (X, Y, Z) 
and non-inertial heliocentric frame of reference (x, y, z) configu¬ 
ration at the present time. 


4.2 From proper motion to the state vector 

We present a method to compute the state vector into the 
inertial Galactocentric frame of reference from proper mo¬ 
tions, radial velocities and distances. 


4-2.1 Absolute proper motion vector 

Since we measure the proper motion of the cluster with re¬ 
spect to the bulge field stars (relative proper motion), in 
order to obtain the absolute proper motion of the cluster 
members we have to correct for the internal proper motion 
of the field stars. Ortolan! et al. (2011) derived for the cluster 
HP 1 an internal proper motion of the bulge field stars equal 
to vt = 1.03 kms“^. Zoccali et al. (2001) assumed to mea¬ 
sure the background field stars, having those in the near side 
moving on average in the opposite direction, with respect to 
those in the far side, because of bulge rotation, making that 
their overall proper motion is null. Terndrup et al. (1998) de¬ 
rived an internal proper motion (/ii, ^b) = (0.2,0.2) mas/yr. 
We concluded that the correction for the internal proper 
motion is small if compared to the uncertainties associated 
to the proper motion of the cluster, so we will ignore it (see 
Ortolan! et al. 2011). 


4-2.2 Initial state vector 

In order to predict the orbit of the cluster, we need the initial 
state vector, i.e. the initial conditions necessary for solving 
the system of differential equations represented by the equa¬ 
tions of motion. The first step is to obtain the value of the 
velocity components of the cluster in the heliocentric frame 
of reference (f7, F, IF). We used a right-handed coordinate 
system, so that U, V, and W are positive in the direction 
of the Galactic center, Galactic rotation and North Galactic 
Pole (NGP) respectively. We refer to Johnson & Soderblom 
(1987) for a more complete explanation of the procedure 
to obtain the heliocentric velocity vector starting from dis¬ 
tance, proper motion and radial velocity. We adopted the 


value of the velocity in the local standard of rest at the 
Sun’s Galactocentric distance equal to ulsr = 239.7 kms“^, 
according to Model III of Irrgang et al. (2013). The velocity 
components of the Sun with respect to the local standard of 
rest are (Schonrich et al. 2010) 

Uq = 11.1 ± 1 km s~^ 

Vq = 12.24 ± 0.2 km s"^ 

IF© = 7.25 ± 0.5 km s~^ 

In©! = 18.04 km s~^ 


Thus, the Sun moves towards the Galactic center and up to¬ 
wards the NGP away from the plane. It also moves around 
the Galactic center faster than it would if it were on a circu¬ 
lar orbit. The proper motion is then measured in the helio¬ 
centric frame of reference, which is not inertial (see Figure 
[HI). According to Terndrup et al. (1998), the perpendicu¬ 
lar components of the velocity vector {y and z direction) 
expressed into the inertial heliocentric frame of reference is 
given by 

Uci.x = demobs + (8) 

where d© is the heliocentric distance of the cluster, fiobs 
the observed proper motion in Galactic coordinates, df is 
the heliocentric distance of the held stars and v©,© are the 
perpendicular components of the velocity vector of the Sun 
with respect to the local standard of rest. Considering that 
the clusters we are studying are projected on the Galactic 
bulge, we assumed that 

df = Rq (9) 

where Rq is the Galactocentric distance of the Sun. Apply¬ 
ing equation [8] we can write 

'A =-([/ + [/©) 

+ ( 10 ) 

The position vector of the cluster in the Galactocentric 
frame (A, Y, Z) is related to the measured position in the 
heliocentric frame of reference (x, y, z) by 



A = Rq — X 

Y =-y (11) 

Z = z . 


5 ANALYSIS OF ERRORS 

In this section we describe the errors affecting our measure¬ 
ments. 


5.1 Astrometric Errors 

Astrometric errors are the combination of different random 
and systematic errors. We identihed several sources of uncer¬ 
tainty, namely the random centering errors, the chromatic 
errors (Ortolan! et al. 2011) and the field distortion errors. 
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We used images of the cluster NGC 6558 as a test case, but 
the results are valid in a more general context. 


Centering Error 

We considered four images of NGC 6558 from NTT+EMMI 
(1993): 

• V band, 10 s 

• V band, 3 min 

• I band, 7 s 

• I band, 5 min 


We selected the same region of the CCD used for the deter¬ 
mination of the proper motion in each image and matched 
the .als tables with the DAOMATCH and DAOMASTER routines. 
We considered the I band images to evaluate the centering 
error, because they are more populated than the images in 
the V band, leading to a better statistics. Furthermore, for 
most of the studied clusters, we used I band images for the 
derivation of the proper motion. We designated the coordi¬ 
nates of the single star in the first image and in the 

second image {x 2 ,y 2 )i, with i = 1 to Nstars (hereafter i will 
be always intended in this range). The mean coordinates of 
the stars are 


I 


Xi = 

Vi = 


(Xl +X2)i 
2 

(yi + y2)i 

2 


( 12 ) 


The standard deviation associated to this mean value is 


(xi - x)} + (X2 - x)l 


{yi - yfi + (2/2 - 2/)i 


1/2 


1/2 


(13) 


In Figure[T3]we show the error on the radial position of stars 
as a function of the instrumental I magnitude 


am = aR -h (Tyf (14) 

We applied the same procedure for the I band images from 
NTT-I-FFOSC (2012) and the results are shown in Figure 

M 

The order of magnitude of the average value of the as¬ 
trometric error is about 0.1 pixel. Taking into account the 
number of about 1000 stars used for the derivation of the 
proper motions of the clusters, the order of magnitude of 
the centering error is as small as 10“® pixel. 


Chromatic Error 

Chromatic errors are due to differential chromatic refrac¬ 
tion (DCR), both atmospheric and instrumental, on stars 
with different colours at different air masses. The refraction 
dependence on wavelength (inverse quadratic to the first 
approximation) makes this effect more pronounced in the 
(V — I) colour than in the near infrared. In order to study 
the displacement of the stars as a function of their colours, 
we subdivided the sample of stars into “red” and “blue” 
groups, evaluating the average astrometric error for both. 
We show the results for NGC 6558 (NTT-I-EMMI 1993) in 
Figure [TS] Note that the number of stars is lower than that 
of the previous figures as a consequence of the match of the 
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Figure 13. Centering error vs. 1 band instrumental magnitude 
for the 1993 images (NTT-I-EMMI). The dashed line shows the 
average astrometric error. 
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Figure 14. Centering error vs. 1 band instrumental magnitude 
for the 2012 images (NTT-I-EFOSC). The dashed line shows the 
average astrometric error. 


I and V .als tables. 

We associate to the chromatic error the quantity 

achr = |dc,r - uc.bl (15) 

We expect to obtain a similar result for the 2012 images. 
It is important to notice that the chromatic effect is small 
if compared with the centring error. In fact, the order of 
magnitude is 10“"^ pixels, and then negligible. 

We performed a further test in order to investigate the 
effects of DCR on our measured proper motions. We se¬ 
lected two clusters of our sample with detected blue HB 
stars, namely NGC 6540 fFigure [23ll and AL 3 (Figure [24l) . 
In order to maximize the colour extension of the selection, 
we selected identified cluster members in a F range that con¬ 
tains both HB and RGB stars (15.5 < V < 16.5 for NGC 
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Figure 15. Centering error vs. instrumental {V — I) colour. The 
red line shows the average centering error for the “red” group, 
while the blue line shows the average centering error for the “blue” 
group; these groups are separated by black vertical line. 
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6540 and 14.5 < < 17 for AL 3). We then analysed the 

value of the proper motion as a function of the star colours. 
Figure 116] shows the results of this analysis. We fitted a re¬ 
gression line (continuous line) to the data points applying a 
linear least square interpolation and plotted the regression 
lines (dashed lines) associated to the 95% confidence bound 
on the coefficient values of the fit model. The right ascen¬ 
sion component of NGC 6540 shows the highest dependence 
of the proper motion on the stars’ colour. The declination 
components of NGC 6540 and both the right ascension and 
declination components of AL 3 show a very low dependence 
on the colour. However, if we consider the uncertainties as¬ 
sociated to the regression model, we note that the all the 
fitted lines are consistent with no slope within the errors. 
We can conclude that we did not find any significant depen¬ 
dence of the proper motion on the star colours and hence 
that we can confidently exclude a strong DGR effect on our 
data. The same conclusion was reached in Ortolani et al. 
( 2011 ). 



Figure 16. Components of the proper motion as function of the 
star colours for NGC 6540 (two two panels) and for AL 3 (bottom 
two panels). 


sian curve and coincides with the difference between the best 
fitting value and the extremes of the 95% confidence bounds. 

Following Johnston et al. (1987), we have all the ingre¬ 
dients to estimate the uncertainties in the velocity compo¬ 
nents. By applying the error-propagation to the estimate of 
the velocity components we found 


5.2 Gaussian Interpolation Error 

The error associated to the double Gaussian interpolation is 

< _ ( 16 ) 

where A 2 ;ci, Aifleid, Ai/d, Ai/fleid are the uncertainties in 
the position of the cluster and field stars centroids in the 
(AX, AY) diagram for the x and y coordinates respectively. 
The uncertainties associated to the astrometric and numer¬ 
ical errors are small if compared with the uncertainties as¬ 
sociated to the Gaussian interpolation, then we can neglect 
them and assume that the main source of uncertainty derives 
from this fitting procedure. This error is estimated directly 
from the routine that computes the coefficients of the Gaus- 





1 


bl2 ■ 6i3 

4 

= C 


~h (/^CK *^77 /tt) j 

I *^7r 

TT^ 

622 • 623 

yw. 



+ (M5Cr,r/7r)^]_ 


532 • b33_ 


(17) 

The elements of the matrix C are the squares of the individ¬ 
ual elements of B, i.e. dj = bfj for all i and j (see Johnston 
et al. 1987). 

The final step to obtain the uncertainties associated to 
the state vector is to propagate the errors in position and 
velocity in the inertial Galactocentric frame of reference. 
This can be easily done computing the partial derivatives of 
the transformations of coordinates described in Section 4. 






































12 Rossi et al. 


6 RESULTS FOR INDIVIDUAL GLOBULAR 

CLUSTERS 

We derived the proper motion of ten bulge clusters. Table [3] 
gives the proper motion values, velocity components of the 
Galactocentric velocity vectors, and other relevant informa¬ 
tion. The input distances and coordinates for the sample 
clusters are from Bica et al. (2006). The data show that the 
proper motions and the velocity components of our targets 
are, on average, lower than those of halo clusters (e.g. Di- 
nescu et al. 1999), with the exception of the two most metal- 
poor clusters (Terzan 4, Terzan 9). For the more reddened 
clusters Terzan 1, Terzan 4 and Palomar 6, we employed the 
z-band (instrumental), which provided a better definition to 
the CMDs. The studied clusters are discussed individually 
as follows. 

Terzan 1 

Terzan 1 is a moderately metal-poor ([Fe/H]~ —1.0) globu¬ 
lar cluster with a dominant red HB. There is evidence that it 
is a second parameter globular cluster showing a post-core 
collapse structure (Ortolan! et al. 1999a). We used as first 
epoch CCD camera on the Danish telescope observations 
taken in June 1990, and as second epoch NTT-1-EFOSC2 
data observed in May 2012, with a baseline of 22 years. 

In Figure[T7]are shown the derived proper motions, the 
double-Gaussian interpolation of the displacements and the 
CMD cleaning of the cluster. The present analysis showed a 
small proper motion both in right ascension and declination, 
that could be an indication of the very elongated orbit of the 
cluster. This picture would be consistent with the relatively 
high radial velocity of Terzan 1 {vr = 114±14 km s”*^, Idiart 
et al. 2002). The values are reported in Table 3. 

Terzan 2 

Terzan 2 is a metal-rich red horizontal branch cluster, lo¬ 
cated very close to the Galactic center (Ortolan! et al. 1997a; 
Valenti et al. 2007). With a baseline of 18 years, a very 
good seeing and a good sampling in the first epoch ob¬ 
servations, taken at NTT-fSUSI, in May 1994, the proper- 
motion-cleaning is very satisfactory. The second epoch was 
taken with NTT-I-EFOSC2 in June 2011. The cluster GMD 
given in Fig. [18] is not richly populated, but the elimination 
of field sequences was effective. Despite its location being 
among the closest objects to the Galactic center, it shows 
a very low proper motion, below 1 mas/yr, together with a 
radial velocity of 109 km/s. 

Terzan 4 

Terzan 4 is among the most metal-poor clusters in the inner 
bulge (Ortolan! et al. 1997b; Origlia et al. 2004). First and 
second epoch observations were obtained respectively, at the 
NTT-bSUSI in May 1994, and NTT-fEFOSG2 in May 2012. 
Both these sets of images have been obtained under excep¬ 
tionally good seeing. The proper-motion-cleaned GMD with 
a baseline of 18 years provided a very well defined RGB. 
The definition of the CMD is comparable to that of NIC- 
MOS/HST one (Ortolan! et al. 2007), showing the possibil¬ 
ity of a few blue HB stars. The proper motion values, as 


well as the Galactic velocities ( lower panel of Table O are 
somewhat higher than those of the sample average, possibly 
characterizing a halo-originated cluster (given also its much 
lower metallicity). The orbit prediction will bring more con¬ 
straints. 

Terzan 9 

Ortolani et al. (1999b) estimated a metallicity of [Fe/H]~ 
—2.0, whereas Valenti et al. (2007) indicates [Fe/H]~ —1.2. 
In the latter case, it could be added to the group of moder¬ 
ate metallicy and BHB. Harris (1996, 2010 version) reports 
a radial velocity of 59 km/s. First and second epoch observa¬ 
tions were obtained respectively, with the DFOSC@Danish 
in July 1998, and and NTT-I-EFOSC2 in May 2012. For this 
cluster the time baseline is of 14 years. It has a relatively 
high proper motion (Table 3). Only a detailed orbital study 
can clarify if it is a halo or a bulge cluster. 

NGC 6522 

NGC 6522 is a moderately metal-poor globular cluster pro¬ 
jected on the Baade’s Window. We derived the proper mo¬ 
tion from two sets of observations: the first epoch images 
have been obtained with the Danish Telescope in June 
1992, and HST/WFPG2 in September 1995, while the sec¬ 
ond epoch images come from NTT (May 2012). Both di¬ 
agrams are well-cleaned, with HST first epoch including 
fainter stars. The dispersion of the proper motion plots is 
tighter for the Danish data, possibly due to the inclusion of 
fainter stars in the HST comparison. The abundance anal¬ 
ysis and radial velocity of 8 stars was presented in Bar- 
buy et al. (2009, 2014). Its proper motion has been inves¬ 
tigated also by Terndrup et al. (1998), who derived values 
with respect to a star at rest in the nuclear bulge equal 
to (/ii, fj.b) = (1-4 ± 0.2, —6.2 ± 0.2) mas yr“^. The top 
6 panels of Figure [21] show the cluster proper motion ob¬ 
tained in the present work. Our analysis produced the val¬ 
ues (/i* , Us) = (3.35 ± 0.60, —1.19 ± 0.34) mas yr“’^, cor¬ 
responding to a proper motion in Galactic coordinates of 
(pi, fib) = (0.4,—3.1) mas yr“^. Our result is about a 
factor of 2 smaller than the one obtained by Terndrup et 
al. (1998). We also tried to derive its proper motion inde¬ 
pendently, by using archive HST images. The results are 
shown in the bottom six panels of Figure !^ This new esti¬ 
mate of the proper motion of the cluster produced the value 
ifia, fis) = (3.6 ± 0.6, —1.2 ± 0.3) mas yr“^, which is con¬ 
sistent with our previous analysis. Note that the calibration 
star from Terndrup is outside our field. We finally also note 
that in Barbuy et al. (2014) a revised radial velocity of Vr = 
-14.3±0.5 km s“^ was obtained, which is basically compati¬ 
ble with the adopted value given in Table [3] 

NGC 6558 

NGG 6558 has a well-determined low proper motion, com¬ 
bined with a high radial velocity of -197.2 km s“^ (Bar¬ 
buy et al. 2007). First and second epoch observations were 
obtained respectively, at the NTT-fEMMI in June 1993, 
and NTT-I-EFOSG2 in May 2012. The second epoch was 
obtained under exceptionally good seeing conditions. The 
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present decontaminated CMD is very populous with well- 
defined sequences. The cleaning was effective in this case, 
superseding previous CMDs for this cluster. Vazques et al. 
(2013) found a value of the proper motion (/r; cos(b),^6) 
= (0.30±0.14, -0.43±0.13) mas/yr. The present proper mo¬ 
tion in Galactic coordinates is (fii cos(b),/r(,) = (0.37±0.43, 
0.29T0.62) mas/yr. Considering the large errors relative to 
the small motion of the cluster, and that the telescopes and 
the techniques are different, the results are in a very good 
agreement. An orbit derivation is needed in order to estab¬ 
lish its dynamics. 


NGC 6540 

NGC 6540 has a detectable small proper motion, and the 
cleaning is effective, showing that we correctly identified 
the CMD features. A very low radial velocity of 17.2 km 
s“^ is given in Harris (1996, 2010 version). First and sec¬ 
ond epoch observations were obtained respectively, with the 
direct CCD camera at the Danish telescope in June 1990, 
and NTT-I-EFOSC2 in May 2012. The observing runs were 
the same as for Terzan 1. Given its low proper motion and 
radial velocity, the object could be near its apogalacticon. 
Bica et al. (1994) estimated a metallicity of [Fe/H]~ —1.0. 
This could be another moderate metallicity cluster with a 
BHB. The proper-motion-cleaning was not very effective, as 
shown by about the same number of BHB stars in the clus¬ 
ter and in the field, due to a very small difference between 
the main motion of the cluster and field stars. Nevertheless, 
the proper-motion dispersion of field stars is sizably larger 
than that of cluster members. We can clearly see that the 
CMD of proper-motion-selected stars (middle panel on the 
bottom of Figure 23) is fairly free of field stars. Most of the 
field stars are indeed included in the bottom-right panel, 
where they should belong. 


AL 3 

First and second epoch observations were obtained respec¬ 
tively, with DFOSC@Danish telescope in March 2000, and 
NTT-I-EFOSC2 in May 2012. AL 3 shows a rather high right 
ascension proper motion as compared with the other cen¬ 
tral clusters, but it is comparable to that of the bulge edge 
cluster NGC 6652. The proper motion in declination on the 
other hand is clearly higher in NGC 6652 than in AL 3. A 
radial velocity measurement for AL 3 is needed for a kine- 
matical diagnostic. As compared with Ortolani et al. (2006), 
a few stars in the RGB are missing in the present CMD, ow¬ 
ing to saturation effects in the second epoch image. The 
proper-motion-cleaned CMD confirms that AL 3 is a BHB 
cluster, and the RGB indicates a moderate metallicity. 


ES0456-SC38 

First and second epoch observations were obtained respec¬ 
tively, with NTT-I-SUSI in May 1994, and NTT-1-EFOSC2 
in May 2012. This cluster has a relevant proper motion. 
Despite this, the proper motion cleaning did not produce 
a well-defined CMD, because of the loose structure of the 
cluster and a high differential reddening. Preliminary radial 


velocity values of ~-160 km s ^ were obtained by S. Vasquez 
(in preparation). 

Palomar 6 

The proper motion of Palomar 6 was based on the longest 
time base of our sample, spanning 20 years. First and sec¬ 
ond epoch observations were obtained respectively, with 
NTT-bEMMI in June 1993, and NTT-fEFOSC2 in May 
2013. The second epoch was obtained with the highest avail¬ 
able spatial resolution. It has a proper motion of around 3 
mas/year. Harris (1996) reports a radial velocity of 190 km 
s“^, and Lee et al. (2004) measured 180 km s“^. The clus¬ 
ter is very reddened with E(B-V)=1.33, and its distance is 
compatible with a location near the Galactic center. For this 
reason, high velocity vectors are expected. The proper mo¬ 
tion cleaning was very effective (see lower panels in Fig. l26ll . 
disentangling the cluster from the bulge population. The 
proper-motion-cleaned CMD confirms that it is a 47 Tuc- 
like cluster, with a red HB (Lee et al. 2004). 

NGC 6652 

We added NGC 6652 to the sample for comparison purposes. 
This borderline bulge cluster is located at 1=1?53, b=-ll?38. 
The edge of the bulge was defined to be at around 10° by 
Soto et al. (2014) and Rojas-Arriagada et al. (2014). Proper 
motions were derived from HST observations as reported by 
Sohn et al. (2014). We find that its kinematics is very differ¬ 
ent from the sample bulge clusters, because its orbit is much 
more extended in the halo, reaching apocentric distances up 
to 10 kpc from the Galactic center. 


7 DISCUSSION AND CONCLUSIONS 

We derived absolute proper motion values for ten inner bulge 
globular clusters, thanks to first epoch images taken by our 
group in the 90s, combined with recent dedicated observa¬ 
tions at the ESO NTT telescope. The proper motions were 
computed from second epoch distortion-corrected positions 
using field stars as absolute reference frame. 

The derived velocity components are as a rule very low 
suggesting that most of these clusters are trapped in the 
bulge and/or the bar. These clusters could have been formed 
earlier than the bar, and were later trapped inside it (Babu- 
siaux et al. 2014). The low velocities are made evident, by a 
comparison with NGC 6652 (Sohn et al. 2014). This cluster 
is located near b~10°, a region considered by Soto et al. 
(2014) and Rojas-Arriagada et al. (2014) as the outer bulge. 
The bulge-halo transition is defined as the region where some 
bulge rotation is still detected for the field stars, but the 
ratio of metal-rich/metal-poor field stars components is in¬ 
verted, with metal-poor stars (halo) prevailing outside this 
region. The comparison of the kinematics between our inner 
bulge sample clusters and NGC 6652 shows a clear differ¬ 
ence in space velocities, indicating the central clusters to be 
a distinct population relative to the outer bulge/inner halo 
clusters. It would be important to acquire proper motions 
of other clusters in the external regions of the bulge. 

Dinescu et al. (1999) found lower velocities for inner 
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galaxy clusters as compared with outer halo ones. We found 
even lower velocities for inner bulge clusters. 

In conclusion, the proper motion analysis of the present 
sample clusters indicates that the moderately metal-poor 
and metal-rich clusters appear to have preferentially low ve¬ 
locity components. The halo contamination inside the bulge 
appears to be limited to very few objects, if any. The proper- 
motion-cleaned CMDs confirm a typical pattern for the Blue 
Horizontal Branch moderate metallicity clusters. 

The orbits and their evolution will be presented in a 
forthcoming paper. 
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CLUSTER 

«(°) 

5(°) 

1 n 

bn 

Ra (mas/yr) 

Mi (mas/yr) 

Vr (km/s) dQ (kpc) 

Terzan 1 

263.95 

-30.47 

357.57 

1.00 

0.51 ±0.31 

-0.93 ±0.29 

114.0 ±14.0 6.2 ±0.6 

Terzan 2 

261.89 

-30.80 

356.32 

2.30 

-0.94 ± 0.30 

0.15 ±0.42 

109.0 ±15.0 8.7 ±0.8 

Terzan 4 

262.66 

-31.59 

356.02 

1.31 

3.50 ±0.69 

0.35 ± 0.58 

-50.0 ±2.9 9.1 ±0.9 

Terzan 9 

270.41 

-26.84 

3.61 

-1.99 

0.00 ±0.38 

-3.07 ±0.49 

59.0 ±10.0 7.7 ±0.7 

NGC 6522 

270.89 

-30.03 

1.02 

-3.93 

3.35 ±0.60 

-1.19 ± 0.34 

-21.1 ±3.4 7.8 ±0.7 

NGC 6558 

272.57 

-31.76 

0.20 

-6.02 

-0.12 ± 0.55 

0.47 ± 0.60 

-197.2 ±1.5 7.4 ±0.7 

NGC 6540 

271.53 

-27.76 

3.29 

-3.31 

0.07 ±0.40 

1.90 ±0.57 

-17.72 ±1.4 3.7 ±0.3 

AL 3 

273.53 

-28.63 

3.36 

-5.27 

4.77 ±0.46 

0.55 ± 0.44 

6.5 ± 0.6 

ESO 456-SC38 

270.45 

-27.83 

2.77 

-2.50 

3.08 ±0.29 

2.00 ± 0.34 

6.7 ±0.6 

Palomar 6 

265.93 

-26.22 

2.10 

1.78 

2.95 ±0.41 

1.24 ±0.19 

181.0 ±2.8 7.3 ±0.7 

NGC 6652 

278.94 

-32.99 

1.53 

-11.38 

4.75 ±0.07 

-4.45 ± 0.10 

-111.7±5.8 9.6±0.9 

CLUSTER 

X (kpc) 

Y (kpc) 

Z (kpc) 

X (km/s) 

Y (km/s) 

Z (km/s) [Fe/H] 

Terzan 1 

1.3 

(2.1) 

0.3 

(0.3) 

0.1 

(0.1) 

-125 ± 14 
(-125) 

-25 ± 18 
(-45) 

-24 ±9 -1.3 

(-24) 

Terzan 2 

-1.2 

(-0.4) 

0.6 

(0.6) 

0.4 

(0.4) 

-122 ± 15 
(-122) 

18 ±22 
(-10) 

-26 ± 14 -0.5 

(-25) 

Terzan 4 

-1.6 

(-0.8) 

0.6 

(0.6) 

0.2 

(0.2) 

46 ± 3 
(46) 

121 ±32 
(92) 

132 ±31 -1.6 

(133) 

Terzan 9 

-0.2 

(0.6) 

-0.5 

(-0.5) 

-0.3 

(-0.3) 

-76 ± 10 
(-76) 

99 ±25 
(54) 

-58 ±16 -1.0 

(-57) 

NGC 6522 

-0.3 

(0.5) 

-0.1 

(-0.1) 

-0.5 

(-0.5) 

20 ±4 
(23) 

-11 ±22 
(-36) 

-129 ± 22 -0.86 

(-127) 

NGC 6558 

0.1 

(0.9) 

0.0 

(0.0) 

-0.8 

(-0.8) 

184 ± 3 
(184) 

-12 ±26 
(-36) 

32 ± 20 -0.97 

(32) 

NGC 6540 

3.8 

(4.6) 

-0.2 

(-0.2) 

-0.2 

(-0.2) 

8±2 

(8) 

-155 ± 19 
(-167) 

20 ±8 -1.0 
(20) 

AL 3 


1.0 

(1.8) 

0.4 

(0.4) 

-0.6 

(-0.6) 

(-) 

(-) 

-1.3 

(-) 

ESO 456- 

SC38 

0.8 

(1.6) 

0.3 

(0.3) 

-0.3 

(-0.3) 

(-) 

(-) 

-0.5 

(-) 

Palomar 6 

0.2 

(1.0) 

-0.3 

(-0.3) 

0.2 

(0.2) 

-195 ±3 
(-195) 

85 ±23 
(62) 

76 ±18 -1.0 

(77) 

NGC 6652 

-1.9 

(-1.1) 

-0.3 

(-0.3) 

-1.9 

(-1.9) 

150 ± 8 
(150) 

169 ± 19 
(138) 

-256 ± 28 -0.81 

(-255) 


Table 3. Kinematical properties. Upper panel: input data; lower panel: derived Galactocentric vectors and velocities, and literature 
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Figure 17. Proper motion and CMD decontamination of Terzan 1. 



Figure 18. Proper motion and CMD decontamination of Terzan 2. 
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Figure 19. Proper motion and CMD decontamination of Terzan 4. 



Figure 20. Proper motion and CMD decontamination of Terzan 9. 
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Figure 21. Proper motion and CMD decontamination of NGC 6522 from Danish Telescope—NTT and NTT-HST data. 
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Figure 22. Proper motion and CMD decontamination of NGC 6558. 



Figure 23. Proper motion and CMD decontamination of NGC 6540. 
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Figure 24. Proper motion and CMD decontamination of AL 3. 



Figure 25. Proper motion and CMD decontamination of ES0456-SC38. 
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Figure 26. Proper motion and CMD decontamination of Palomar 6. 













